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ABSTRACT 

JHKs photometry is presented from a three-year survey of the central regions of 
the Local Group dwarf irregular galaxy IC 1613. The morphologies of the colour- 
magnitude and colour-colour diagrams are discussed with particular reference to the 
supergiants and M- and C-type asymptotic giant branch (AGB) stars. Mean JHKs 
magnitudes, amplitudes and periods are given for five O-rich and nine C-rich Mira 
variables for which bolometric magnitudes are also estimated. A distance of 750 kpc 
((m — M) o = 24.37 ± 0.08 mag) is derived for IC 1613 by fitting a period-luminosity 
relation to the C-rich Miras. This is in agreement with values from the literature. 
The AGB stars exhibit a range of ages. A comparison with theoretical isochrones 
suggests that four luminous O-rich Miras are as young as 2 x 10 s yrs. One of these 
has a lithium absorption line in its spectrum, demonstrating that it is undergoing 
hot bottom burning (HBB). This supports the idea that HBB is the cause of the 
high luminosity of these AGB stars, which puts them above the fundamental period- 
luminosity (PL) relation. Further studies of similar stars, selected from their positions 
in the PL diagram, could provide insight into HBB. A much fainter, presumed O- 
rich, Mira is similar to those found in Galactic globular clusters. The C Miras are 
of intermediate age. The O-rich variables are not all recognized as 0-ricli, or even as 
AGB stars, on the basis of their J — K$ colour. It is important to appreciate this when 
using near-infrared surveys to classify AGB stars in more distant galaxies. 

Key words: stars: variables: AGB; stars: carbon; galaxies: distances and redshifts; 
galaxies: individual: IC 1613; (galaxies:) Local Group; infrared: stars 


1 INTRODUCTION 

IC 1613 is a Local Group dwarf irregular galaxy. Its prox¬ 
imity (0.72 Mpc Scowcroft et al. 2013) and isolation, to¬ 
gether with its low foreground extinction (e.g. foreground 
E(B — V) = 0.025 mag; Schlegel et al. 1998) have made it an 
attractive target for various programmes. Studies of young 
stars show a metallicity similar to that of the SMC, but with 
a low oxygen abundance (Garcia et al. 2014; Bouret et al. 
2015), while older stars appear to be more metal-deficient 
(e.g. Skillman et al. 2014). A recent investigation of the star 
formation history (Skillman et al. 2014) supports earlier sug¬ 
gestions that it has been constant over time, without any 
dominant phase of star formation. 

Borissova et al. (2000) studied the distribution of lu¬ 
minous cool stars from J and A'-band imaging and found 
AGB stars covering a wide range in age. Albert, Demers & 
Kunkel (2000) conducted a narrow-band wide field survey 
of IC 1613 for C and M stars, and found the (intermediate 
age) C stars extended out to 15 arcmin, well beyond the 


regions where star formation is currently active. Bernard et 
al. (2007) conducted a wide field optical survey of IC 1613, 
and traced red giant branch (RGB) stars out to a radius of 
16.5 arcmin (3.6 kpc), showing the galaxy to be more ex¬ 
tended than previously thought. Sibbons et al. (2015) use 
JHKs photometry to identify the tip of the RGB (TRGB) 
at Kq = 18.25 ± 0.15 mag and to discuss the relative num¬ 
bers of O- and C-rich AGB stars. Jackson et al. (2007) and 
Boyer et al. (2009; 2015) use Spitzer mid-infrared photom¬ 
etry to discuss dust production and demonstrate that there 
are AGB stars in IC 1613 that are too red to be found by 
surveys at shorter wavelengths. 

Various supergiants, including Cepheids and other vari¬ 
ables, were identified by Sandage (1971). Spectra were taken 
of some by Humphreys (1980) and infrared photometry ob¬ 
tained by Elias & Frogel (1985). More recently surveys for 
variables have been made by Saha et al. (1992), Antonello 
et al. (1999, 2000), Mantegazza et al. (2001), Bernard et al. 
(2007, 2010) and OGLE II (Soszyriski et al. 2009). 

The present work forms part of a broad study of AGB 
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variables in Local Group Galaxies which so far has cov¬ 
ered the dwarf spheroidals: Leo I (Menzies et al. 2002; 
2010), Phoenix (Menzies et al. 2008), Fornax (Whitelock 
et al. 2009) and Sculptor (Menzies et al. 2011), as well as 
NGC6822 (Whitelock et al. 2013), another dwarf irregular. 


2 OBSERVATIONS 

We used the Japanese-South African IRSF equipped with 
the SIRIUS camera, which permits simultaneous imaging 
in the J, H and Ks bands (see Nagayama et al. (2003) 
for details). We defined 4 overlapping fields, each approx¬ 
imately 7.8 arcmin square, with field 1 centred at a(2000.0) 
= 01:05:01.3 and <5(2000.0) = 02:11:43.8. Field 2 is 7 arcmin 
W of field 1, fields 3 and 4 are 7 arcmin south of fields 2 and 
1, respectively. Fields 1 and 2 were observed in JHKs at 
11 epochs, and the other two at 10 epochs, over a period of 

3 years. Typically 30 dithered exposures of 30 s each were 
combined at each epoch, though occasionally, depending on 
sky brightness at Ks, exposure times were reduced to 20 s. 

The basic data for stars with standard errors less than 
0.1 mag in each band are provided on-line, and the first few 
lines of the catalogue are illustrated in Table [I] (The Mira 
variables, discussed in section 8, are not in this table). 

Observations were obtained with the Robert Stobie 
Spectrograph on the Southern African Large Telescope 
(SALT). An 1800 line/mm grating was used to cover the 
wavelength region from 6000 to 7200 A with a resolving 
power of about 4000. Two spectra with exposure times 
of 2400 s each were obtained at GJD 2456980.3903 and 
2456981.3737, respectively and reduced with the SALT 
pipeline (Crawford et al. 2010). 

J?-band photometry was obtained with the Las Cumbres 
Observatory Global Telescope (LCOGT) Network. This was 
primarily for the purpose of establishing when the M-type 
Mira, G3011, was sufficiently bright to obtain spectroscopy 
with SALT (see section 4), but it also served to give us 
some variability information on all four of the bright M-type 
Miras. 

The reddening to IC 1613 is small and no corrections are 
applied for the various diagrams. For the bolometric mag¬ 
nitudes a small correction is applied as discussed in section 

9. 

A mosaiced Ks image is shown in Fig. [l] 


3 UKIRT DATA 

Infrared photometry of IC 1613 has been obtained with 
UKIRT and published by Lawrence et al. (2007), as part 
of the UKIDSS survey, and by Sibbons et al. (2015). It is 
potentially useful to combine these data with our IRSF mea¬ 
surements in order to determine periods for the large ampli¬ 
tude variables. 

The Sibbons et al. (2015) photometry (Jw, Hw, Kw) 
was re-reddened assuming E ^ B _ V ) = 0.02 so that 
Aj,Ah,Ak are 0.015, 0.009, 0.005 mag, respectively, then 
corrected to the Two Micron All Sky Survey (2MASS) sys¬ 
tem {J 2 ,H 2 ,K 2 ) using: 

(1) 


H 2 = Hw + 0.081(Jw — Hw) — 0.033 (2) 

K 2 — Kw — 0.010(Jw — Kw) (3) 

which are derived from Hodgkin et al. (2009). 

The UKIDSS data from Lawrence et al. had not been 
reddening corrected and were therefore simply transformed 
using the above equations. 

As discussed in our paper on NGC6822 (Whitelock et 
al. 2013) IRSF and UKIRT photometry for non-variables 
shows good agreement for H and Ks, but a colour term for 
J of the order of 0.1 mag for stars with J — Ks = 1.0, in the 
sense that the transformed Sibbons et al. J mag is fainter 
than our values. 

We assume that the above transformation applies to the 
large amplitude variables, while noting that this may not be 
the case because of their unusual energy distributions. 


4 SPECTRUM OF G3011, THE LUMINOUS M 

GIANT 

The SALT spectrum of G3011 is shown in Fig. [2] together 
with one of HV11329, a super-lithium-rich star in the SMC 
(Smith, Plez & Lambert 1995), taken with the same set up, 
for comparison. The region around the lithium line is shown 
at higher magnification. 

The presence of the lithium doublet at 6708A is quite 
clear as is the ZrO band head at 6474A As an indication of 
the strength of the lithium line, we find an equivalent width 
of about 350 m A, using a pseudo-continuum defined by the 
fluxes at 6704 and 6712A. With the same procedure, the 
equivalent width in our HV11329 spectrum is about 480mA 
(compared with 767mA given by Smith et al. (1995), using 
a ’local flux maximum’ as continuum). The presence of s- 
process elements leads us to classify this as an AGB star 
while the lithium clearly indicates hot bottom processing. 

The heliocentric radial velocity of G3011 based on the 
position of Ho? is —233 Kms -1 , which is very similar to that 
of other stars in IC 1613 (Kirby et al. 2014) as well as to the 
HI velocity (Hoffman et al. 1996). What is perhaps quite 
surprising is the strength of the metal lines in the same 
region, note in particular the line just to the red of lithium 
with a rest wavelength of 6710.5 A which must be mostly 
iron (Davis 1947), and which is not seen at all in the SMC 
star spectrum. This might reasonably lead us to question 
the assumption that IC 1613 is a very metal weak galaxy. In 
that context it is interesting to note that Garcia et al. (2014) 
in looking at the spectra of young hot stars recorded “An 
unexpected finding is that when comparing analogue stars 
in the Fe V and Fe IV pseudo-continua, the line strength of 
the Fe forest in IC 1613 stars is similar to that in SMC stars, 
or even stronger. We expected a significant difference given 
IC 1613’s poorer metal content. ” 

The metallicity of IC 1613 requires a proper abundance 
analysis of stars with atmospheres that can be properly mod¬ 
elled. 


5 COLOUR-MAGNITUDE DIAGRAM 

The colour-magnitude and two-colour diagrams for stars 
from Table [l] plus the variables discussed in section 8, are 


J 2 = Jw + 0.075(JV — Hw) + 0.002 
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Table 1. Data for stars with standard errors less than 0.1 mag (note that this selection omits the large amplitude variables which are 
listed in Tables 4, 6 and 7 below and G1027, the Cepheid variable). The full table is available on-line. The first two columns are the 
equatorial coordinates in degrees; G is our own identification number; the mean photometry, JHKg is listed together with its standard 
deviation, SJHKg. 


RA Dec 

(equinox 2000) 

G 

J 

SJ 

H 

SH 

K S 

(mag) 

SK s 

J - H 

H-K S 

J-K S 

16.21004 

2.17447 

1001 

12.512 

0.025 

12.171 

0.016 

12.140 

0.032 

0.341 

0.031 

0.372 

16.31182 

2.23713 

1002 

13.536 

0.018 

12.954 

0.032 

12.729 

0.039 

0.582 

0.225 

0.807 

16.25696 

2.14422 

1003 

13.934 

0.044 

13.234 

0.045 

13.013 

0.056 

0.699 

0.221 

0.921 

16.24328 

2.15233 

1004 

14.096 

0.032 

13.427 

0.036 

13.191 

0.034 

0.669 

0.236 

0.905 



Figure 1. Mosaic in the Ks band of the fields observed in IC 1613, with north up, east to the left. The positions of various categories of 
stars referred to later in the text are shown as squares of different colour: Supergiants red, O-rich Miras blue, Large amplitude C Miras 
cyan, and other large amplitude variables magenta. 


illustrated in Fig. [3] and Fig. [4] Sibbons et al. (2015) found 
the TRGB to be at Kg = 18.25 ± 0.15 mag. Our limiting 
magnitude is Kg ~ 18.0, so all of the objects in our sample 
that are in IC 1613 (i.e. not foreground stars or background 
galaxies) are either supergiants or AGB stars. 

Following the same procedure as for NGC 6822 (White- 
lock et al. 2013) we assume that stars with J — H < 0.75 
are dwarfs. Note that Sibbons et al. used a somewhat 
bluer colour to make the same division. The identification 
of dwarfs is helpful for an interpretation of the colour- 
magnitude diagram, but it is not a factor in our identi¬ 
fication of AGB variables. However, both the supergiants 


(section 7) and the O-rich AGB variables (section 8) have 
similar colours to dwarfs, so any we have not specifically 
identified may be misrepresented as dwarfs in the figures. 

The colour-magnitude and two-colour diagrams (Figs. [3] 
and [4| are very similar to those of NGC 6822 (figs. 1 and 
2 Whitelock et al. 2013), when we take into account the 
observational limits and the different distances to the two 
galaxies. 

In IC 1613 Fig. |f] there are rather more points with 
0.5 < H — Ks < 1.0 and 0.4 < J — H < 0.8 than 
in NGC 6822. The faint (presumed M-type) AGB variable 
(G4237) falls in this area as do two stars that are tenta- 
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Figure 2. (a) Spectrum of G3011 (top) compared with that of HV11329 (bottom). Both have been scaled to 1.0 at 6600.4, and the 
SMC star’s spectrum is shifted downwards by 0.5 units for display purposes. Ha is strongly in emission in G3011. (b) The region around 
the lithium doublet (marked as a vertical line). 


Table 2. Photometry of objects identified as supergiants in the literature 


G 

RA 

equinox 

Dec 

2000 

J 

H 

K S 

<jJ aH 

(mag) 

crKs 

R 

1027 

16.25295 

2.18007 

16.03 

15.63 

15.52 

0.70 

0.77 

0.78 


4013 

16.22231 

2.08931 

15.464 

14.757 

14.597 

0.052 

0.072 

0.088 

17.551 ± 0.041 

1004 

16.24328 

2.15233 

14.096 

13.427 

13.191 

0.032 

0.036 

0.034 

16.430 ±0.100 

1003 

16.25696 

2.14422 

13.934 

13.234 

13.013 

0.044 

0.045 

0.056 

15.881 ± 0.021 

1009 

16.27144 

2.19849 

15.140 

14.479 

14.315 

0.023 

0.022 

0.035 

16.288 ±0.027 

1010 

16.29004 

2.20715 

14.900 

14.175 

13.985 

0.045 

0.039 

0.040 

16.889 ±0.019 

1008 

16.31963 

2.18765 

14.834 

14.117 

13.923 

0.037 

0.056 

0.036 

17.000 ±0.024 

3003 

16.16047 

2.01605 

14.161 

13.512 

13.342 

0.025 

0.086 

0.029 

16.303 ±0.014 

3188 

16.16983 

2.01953 

18.16 

17.73 

16.78 

0.16 

0.09 

0.08 

19.234 ±0.022 

1025 

16.25852 

2.17354 

16.693 

15.944 

15.702 

0.084 

0.062 

0.066 

18.422 ± 0.055 


The a values for the Cepheid, G1027, are the peak-to-peak amplitudes. 


tively identified as C-rich (G2143, G2075) from the narrow- 
band photometry (see section 6). It is certainly possible that 
these particular stars suffer some weak contamination in the 
J band. In Fig. [ 3 ] most of the rest of the points fall among 
the bluer C stars, but have not been identified as such. It 
is likely that these sources are unresolved distant galaxies. 
Note also that G3188, discussed in section 7, has the most 
extreme colours of this group and is probably a background 
source. 


6 DISTINGUISHING C- AND M-TYPE AGB 
STARS 

Narrow band CN and TiO photometry for IC 1613 by Albert 
et al. (2000) was used to identify M and C stars in our 
sample. We considered only stars for which the uncertainty 


on R— I is a(R — I) < 0.09 mag and used the colour criteria 
of Albert et al. (see their fig. 5). These stars are identified 
in the colour-magnitude (Fig. [3| and two-colour (Fig. [4| 
diagrams. 

Sibbons et al. examined stars brighter than the TRGB. 
They divided their sample into C and M stars on the basis of 
the J — Ks colour, such that stars with 2MASS (J — Ks) < 
1.22 mag are probably M stars and those with 2MASS (J — 
Ks) > 1.22 mag, are probably C stars. They also noted 
that there was some overlap at the dividing line. In Fig. [5] 
we show a colour-magnitude diagram based on photometry 
from Sibbons et al., in the area surveyed by us (except for the 
variable stars which are IRSF measurements). The Sibbons 
photometry has been put on the 2MASS system (section 3), 
so this diagram can be compared directly with Fig. [3] The 
two diagrams are similar except for the fainter limits of the 
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Table 3. Other identifications for the supergiants 


G 

Sandage 

OGLE 

V V-I 

(mag) 

GCVS 

Sp 

Sibbons 

Other 

Ref 

1027 

V22 




V0045 DCEP 

K2I 


P=146.35 


4013 

V23 

11437 

18.359 

1.844 

V0038 Lc 


10977 

V2946 P=645 

1 

1004 

V32 

13677 

17.409 

2.056 

V0042 Lc 

Mia 

15855 


2 

1003 

V38 

13675 

17.324 

1.987 

V0046 Lc 

MOIa 

15299 


2 

1009 

V43 

18839 

17.662 

1.538 

V0053 Lc 

Mil 

18518 


2 

1010 

V58 

18838 

17.666 

1.677 

V0056 Lc 


18909 



1008 

V45 

17879 

17.672 

1.740 

V0060 SRc 


17962 

P=529 

GCVS 

3003 






M2-4I 

7542 


3 

3188 


3692 

19.326 

0.186 


M0-2I 


V0037D 

3, 4, 5 

1025 

V39 

13682 

18.851 

1.234 


sg B[e]/LBV? 

17209 

OGLE P=14.33 

6 


Ref: 1 period Antonello et al. (2000); 2 Sp type Humphreys (1980); 3 Britavsky et al. (2014); 4 Mantegazza et al. (2001); 

5 see text; 6 Herrero et al. (2010) 



Figure 6. Light curve for the Cepheid variable G1027 with each 
point plotted twice to emphasize the variability. The IRSF and 
UKIRT observations are illustrated as closed and open circles 
respectively. The line is a third order sine curve with P=145 days. 

UKIRT data, so they see stars which are both fainter and 
redder than we do. There is no major systematic difference 
between the Sibbons et al. and the Albert et al. divisions 
into C and M stars, and the Albert et al. results indeed 
indicate that there is no strict separation at J — K = 1.2, 
but a good deal of overlap. 


7 SUPERGIANTS 

It is important to recognize the M supergiants in IC 1613, 
particularly the variables, for the purpose of comparison 
with the luminous AGB stars discussed below. 

We also clearly identify the long period Cepheid 
(G1027) which was discovered by Sandage (1971), as V22, 
with a period of 146.35 days, and who noted that it did not 


repeat well. Antonello et al. (1999) call the same star V2396 
and find P=145.6 days, but do not have enough data to cover 
the full cycle. Our period of 147 days from the IRSF data 
alone, or 145 days including the two UKIRT observations 
(see Fig. [6|, is consistent with both these values. The same 
star (V22) in Scowcroft et al. (2013) was assigned a period 
of 123.88 days taken from an earlier paper where they had 
insufficient data to determine the period. This period does 
not fit any of the available data. There are of course many 
more, shorter period Cepheids in IC 1613 (e.g. Udalski et al. 
2001 - henceforth OGLE II; Scowcroft et al. 2013) which are 
too faint to show up in this survey. 

In addition to the bright Cepheid, Table [2] contains 
mean values of the IRSF JHKs and LCOGT R-band pho¬ 
tometry of the stars which are identified as supergiants in 
the literature (see Table [3|. The mean photometry is based 
on at least 10 observations at JHKs and 15 LCOGT ob¬ 
servations for the R magnitudes. For these bright stars the 
standard deviation, cr, gives an indication of the level of 
variability, which is small in all cases (other than for the 
Cepheid). It is certainly too small to look for periodicities, 
although most of these stars have been identified as vari¬ 
ables, e.g. by Sandage (1971), as indicated in the second 
column of Table [3] The OGLE II identity, magnitude and 
colour are given in columns 3, 4 and 5, respectively; the 
General Catalogue of Variable Stars (GCVS - Samus et al. 
2013) identity and variability class is given in column 6, the 
spectral type in column 7 and the Sibbons et al. catalogue 
number in column 8. Column 9 contains any other variabil¬ 
ity information while column 10 contains the reference to 
this information and/or to the spectral type. 

With the exception of the Cepheid, G1027, and G3188 
the supergiants have very similar colours and in Fig. [I] ap¬ 
pear on the border between the M giants and the dwarfs, 
but distinct from the luminous M giant variables discussed 
below. 

Herrero et al. (2010) investigated G1025, which they 
describe as a Luminous Blue Variable (LBV). Its colours in 
pig. a are similar to the late-type supergiants, but is even 
fainter in Fig. [3] than the bright M-type variables. It is pos¬ 
sibly a binary or a confused image. 

G3188 is included in the tables because it was identi¬ 
fied as IC1613-1 by Britavsky et al. (2014) with a spectral 
type of M0-2I. However, Mantegazza et al. (2001) identify 
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Figure 3. Colour-magnitude diagram showing the stars from the various tables. Red and blue (or cyan) symbols represent stars with 
photometry from Albert et al. (2000) that shows them to be C or M stars, respectively. Variables are shown as green circled in black. 
Stars with J — H < 0.65 are probably dwarfs and are shown in gray, or cyan if the Albert et al. photometry suggests they are M stars. 
Those Miras for which we can estimate periods are shown as asterisks, while other large amplitude variables are illustrated as triangles. 
Supergiants are large magenta closed circles. Black symbols are used where there is no information about spectral type, no indication 
of large amplitude variability and no reason to think the star is a dwarf. The lines represent isochrones from Marigo et al. (2008) at 
a distance modulus of 24.37 mag; the solid line and dotted lines represents populations of age 10 9 yrs with metallicity Z=0.001, and 
age 2.5 X 10 s yrs with Z=0.004, respectively. Stars on the blue and red parts of the tracks are expected to have C/O < 1 and > 1, 
respectively. 


the same object as a blue variable, V0037D, with a period 
of 62.71 days. The OGLE II data confirm its blue colour, as 
does the Sloan Digital Sky Survey (r — i = —0.16 mag; Ahn 
et al. 2012). Our infrared images show it as extended and 
the UKIDSS catalogue (Lawrence et al. 2007) indicates that 
it has a high probability of being a galaxy. The UKIDSS 
infrared magnitudes are similar to ours. Inspection of the 


Spitzer catalogs and images from Boyer et al. (2015) show 
multiple sources within a three arcsec radius at these coor¬ 
dinates. In Fig. HI it appears among the C stars. On the basis 
of its radio and X-ray emission Flesch (2010) finds there is a 
97 percent chance that the object is a QSO. He gives 0 and 1 
percent chances of its being a galaxy or a star, respectively. 
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H —K s 


Figure 4. Two-colour diagram with the same stars as in Fig. [3] The line represents the locus for Galactic carbon Miras (equation 2 
from Whitelock et al. (2006) converted to 2MASS as: (H — Ks )0 = 1.003(J — H )0 — 0.428. The symbols are the same as in Fig. [3] 


It is not clear what this is, but it is certainly not an M-type 
supergiant and it is probably more than one source! 


8 VARIABLES 

In the following, in addition to our IRSF photometry, we 
make use of the infrared photometry from UKIRT described 
by Lawrence et al. (2007) and by Sibbons et al. (2015), 
converting it to the 2MASS system as described in Sec¬ 
tion 3. These observations were taken between three and 
seven years after our measurements. This cadence is not 
ideal and long term changes in the mean magnitude, which 


are common among C-rich stars, can make determining pe¬ 
riods more difficult. There also remain possibilities of dif¬ 
ferences in the photometric systems for stars with extreme 
colours. Nevertheless, as the light curves illustrated later 
demonstrate, the extra observations do help with many of 
the period determinations. 

The bright Cepheid variable was discussed in section 7 
and its photometry is listed in Table [2] 

We examined the light curves of all stars with Ks < 17 
mag for which we had at least 10 observations and which 
showed a standard deviation in J, H or Ks of > 0.2 mag, 
going to lower standard deviations for brighter magnitudes. 
Given that our primary objective was to find large amplitude 
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Figure 5. Colour-magnitude diagram showing the variables (IRSF observations) together with photometry from Sibbons et al. (2015) for 
the other stars, which has been reddened and converted to the 2MASS system as described in the text. Red and blue symbols represent 
stars that Sibbons et al. classify as (probable) C- and O-rich AGB stars, respectively. 


AGB variables we also examined stars with J — Ks > 2.2 
mag, finding them all to be variable at some level. 

Variables with J—Ks >1.2 mag are defined as probable 
C-type stars, while the bluer ones are assumed to be M-type, 
i.e. O-rich stars (see section 6). Many of these M- and C- 
type classifications are confirmed from their narrow-band 
colours (Albert et al. 2001) or from spectroscopy in the case 
of G3011. 

Where possible, periods were determined by Fourier 
analysis of the H and/or Ks data. For the large amplitude, 
Mira variables for which periods could be determined the 
Fourier mean JHKs magnitudes and peak-to-peak ampli¬ 


tudes, A J, AH, A Ks were also measured. These are listed 
in Tables [4] and [6] for the probable O-rich and C-rich stars, 
respectively. We followed the same practice as in our ear¬ 
lier papers (e.g. Whitelock et al. 2006) of classifying any 
star with a measurable period and an amplitude at Ks of 
SKs > 0.4 mag as a Mira, and anything with a smaller am¬ 
plitude as a semi-regular (SR) variable. The Sibbons et al. 
identification is given and where they identify the source as 
an AGB star a type of M or C is also given. If no M or C is 
given then Sibbons et al. did not recognize the star as be¬ 
ing on the AGB. The ID and type from Albert et al. is also 
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Figure 8. H Light curve for the variable G4200. Closed and 
open circles are IRSF and UKIRT photometry, respectively. 


listed where the uncertainty of their R — I was or-i < 0.09 
mag. 

Note that G3083 is listed twice with two possible peri¬ 
ods, 147 and 280 days, as it was impossible to distinguish 
between these with the available JHKs data. The one year 
alias of 280 days is 158 days, i.e. not exactly 147, and possi¬ 
bly this is a case of aliasing in rather noisy data. Spitzer data 
(Whitelock et al. in preparation) favours a longer period and 
we therefore adopt 280 days in the following discussion. 

Table [ 7 ] contains the photometry for the stars for which 
we could not determine periods. Some of these will be Miras 
with erratic variations or insufficient data, others will be SR 
variables. Most of them are C-rich as determined by Albert 
et al. (2001). The tabulated A Ks provides the approximate 
range of the variability. No attempt is made to list stars with 
amplitudes AK < 0.4 mag, of which there are many. 

Figure [ 7 ] shows the Ks light curves for the stars in Ta¬ 
ble [6] together with their best fit first-order sinusoid at the 
period noted on the figure. The UKIRT observations (open 
circles) are in reasonably good agreement with our data, 
given that they were obtained at a different time and con¬ 
verted to the 2MASS system as described in section 3. This 
is important because Whitelock et al. (2013) noted signifi¬ 
cant differences between their mean Ks mags and those by 
Battinelli & Demers (2011) for the stars in common (see 
their fig. 6), and there was a possibility that this was due to 
transformation problems with the colours of these very red 
stars. 



phase 


Figure 9. Ks Light curve for the faintest presumed M-type 
Mira. Symbols as for Fig. [7] 

In the colour-magnitude (Fig.[3| and two-colour (Fig.[4| 
diagrams the C-rich Miras are seen to be the reddest of the 
stars we observe and we would anticipate that they have 
large mass-loss rates, as the Spitzer photometry (Boyer et 
al. 2015) also suggests. The other variables overlap in colour 
with the Miras, supporting our expectation that some of 
them are Miras. The light curve of G4200, which is the red¬ 
dest of these, is shown in Fig. [8] There is no sign of period¬ 
icity, although it is highly variable over a timescale of about 
3000 days. 

The luminous M variables are discussed separately be¬ 
low. G4237 (at Ks = 17.28 and J — Ks = 1.19 mag) is 
significantly fainter than these objects; it can be seen about 
2 mag below them in the colour-magnitude diagram (Fig. [3] 
Its light curve is shown in Fig. [9] It was not recognized as an 
AGB star by Sibbons et al. because they would have classi¬ 
fied it as a dwarf from its (J — H )o colour. Our mean value 
was J — H = 0.52 mag, so we too would class it as a dwarf 
on the basis of colour. It is difficult to be certain without a 
spectrum, but the light curve and luminosity are most sim¬ 
ply explained if it is a Mira and it is discussed further in 
section 9. 

It is notable that the four bright Miras have similar 
luminosities and periods to each other, and are therefore 
of very similar age. It is therefore significant that there are 
no Miras, and rather few AGB stars, in the two magnitude 
interval fainter than these stars and none brighter. These 
four stars are discussed further below (section 9). 

8.1 Completeness of survey for large-amplitude 
variables 

As in our analysis of NGC 6822 (Feast et al. 2012; Whitelock 
et al. 2013) we use a reference frame in the H band to pro¬ 
vide positions at which “fixed-position” DoPHOT photome¬ 
try was performed. It is therefore possible that extremely red 
stars or stars with a very large amplitude of variation might 
have been missed. If the star was not measurable on the ref¬ 
erence H frame, then it would not have been retained in the 
catalogue even if found in any other band or on any other 
frame. The limiting magnitudes at J, H and Ks for this sur¬ 
vey are approximately 20.3, 18.3 and 18.0 mag, respectively. 
This means that we would have missed red variables with 
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Figure 7. K$ Light curve for the C-type Miras arbitrarily phased (zero at JD2450000); each point is plotted twice to emphasize the 
variability. The best-fitting first order curve is also illustrated. It is from this that the mean magnitude and amplitude were determined. 
Closed and open circles are IRSF and UKIRT photometry, respectively. 
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Figure 10. Kg Light curve for the four luminous M-type Miras; 
symbols as for Fig. [7] 

H — Ks = 2.0 mag that were fainter than about Ks = 16.3 
mag or J = 18.4 mag at the time the reference frame was 
obtained, even though these latter values are significantly 
above the relevant limiting magnitudes. 

Furthermore, it is important to recognize that we do 
not have as good temporal coverage or cadence as we did for 
NGC6822 (Whitelock et al. 2013), IC 1613 is more distant 
than NGC6822 (by between 0.7 and 0.8 mag), and it is 
therefore clear that we do not expect our sample of large 
amplitude variables in IC 1613 to be complete. 

Table [8] lists the stars from Sibbons et al. (2015) with 
J — Ks > 2 mag that are not in our catalogue. All except 
the three faintest ones (20671, 8310 and 20772) are visible on 
our H reference frames, but too faint to measure. Several of 
them have UKIDSS measurements that are a few tenths of 
a magnitude different from the Sibbons et al. values. All but 
one of them, 16335, have Spitzer measurements from Boyer 
et al. (2015) and 16 of those (all except 8310) have J— [3.6] > 

3.1 mag, making them extreme AGB stars according to the 
definition of Blum at al. (2006) and therefore very good 
candidates for large amplitude C-rich variables. 

8.2 Luminous O-rich variables 

The first four presumed O-rich AGB stars listed in Table [4] 
are well separated from the other variables in the colour- 
magnitude diagram (Fig.[3|. Their Ks light curves are illus¬ 
trated in Fig. 1 10t while Fig. |ll| shows their LCOGT i?-band 



0 0.5 1 1.5 3 0.5 1 1.5 3 

phase 


Figure 11. R Light curve for the four luminous M-type Miras 
from LCOGT data. 


observations, phased on the periods determined from the in¬ 
frared photometry. These periods fit reasonably well and in 
all cases show that the peak at R occurs earlier than at Ks, 
as might be expected. The R amplitudes will be significantly 
larger than the Ks amplitudes. Soszyhski et al. (2009) use 
an I amplitude of Al > 0.8 mag to define Miras within the 
very extensive OGLE data and we would expect the R am¬ 
plitude to be slightly greater than that at I. Thus all of these 
stars would be classed as Miras from their R photometry. 
Our knowledge of the individual stars is summarized briefly 
in Table [3] and the details are discussed below. 

The data summarized in Table [3] show that two of these 
stars, G1017 and G1016, are O-rich from their narrow-band 
colours, while a third, G3011, is O-rich from its published 
spectral type, M3III (see also section 4). All four are large 
amplitude variables with J — Kg < 1.2 mag and we therefore 
assume that all four are O-rich Miras. Sibbons et al., work¬ 
ing only with JHKs colours and no variability information, 
found only one O-rich AGB star, G2035, two of the others 
were not classed as AGB stars at all while the fourth, G1016, 
was classed as a C-rich AGB star. Because JHKs photom¬ 
etry has become an important method for identifying AGB 
stars in other galaxies it is vital to appreciate that it does 
not work for these important luminous AGB variables (nor 
even for the fainter presumed O-rich Mira G4237 discussed 
above). Some notes on three of the variables follow; 

G1017: There is a variable in the Catalina survey 
(Drake et al. 2014), J010506.6+021333, that is 2.7 arcsec 
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Table 4. Large amplitude variables (O-rich and probable O-rich) 


G 

RA 

equinox 

Dec 

2000 

P 

(d) 

J 

H 

K S 

J ~ Kg 

(mag) 

A J 

AH 

A K S 

^bol 

1017 

16.27831 

2.22648 

580 

16.25 

15.61 

15.14 

1.11 

0.89 

1.04 

0.95 

18.332 

3011 

16.20428 

2.09140 

550 

16.18 

15.57 

15.17 

1.01 

0.64 

0.90 

0.82 

18.261 

2035 

16.11180 

2.17299 

530 

16.33 

15.67 

15.20 

1.13 

0.75 

0.74 

0.75 

18.402 

1016 

16.26576 

2.19903 

464 

16.18 

15.32 

15.04 

1.14 

0.43 

0.49 

0.50 

18.085 

4237 

16.23002 

2.07187 

178 

18.47 

17.95 

17.28 

1.19 

0.79 

0.74 

1.07 

20.595 


Table 5. O-rich variables - other data 


G 

Sibbons+ 

Albert+ 

OGLE II 

/ V-I 

(mag) 

SP/Var 

ref 

1017 

19757 - 

3391 M 

18862 

19.584 

2.989 



3011 

11122 - 


4725 

18.932 


M3III 

Kurtev+2001 

2035 

17180 M 


1966 

17.729 

2.887 

VC208 

Bernard+2010 

1016 

18536 C 

529 M 

15569 

17.312 

2.292 

V40,V1872,V51 

1 ,2, GCVS 

4237 

9825 - 


10360 

19.942 





Ref: 1 Sandage (1971), 2 Antonello et al. (1999). 


Table 6. Large amplitude variables (C-rich and probable C-rich) 


G 

RA 

equinox 

Dec 

2000 

P 

(day) 

J 

H 

K S 

J-K S 

(mag) 

A J 

AH 

A K S 

T^bol 

Sibbons+ 

Albert+ 

1147 

16.28988 

2.23202 

879 

18.54 

17.06 

15.64 

2.90 

1.34 

1.36 

1.04 

18.729 

19976 C 


3235 

16.18318 

2.09688 

425 

18.92 

17.34 

16.12 

2.80 

0.92 

0.78 

0.52 

19.363 

11555 C 


4183 

16.20909 

2.11792 

410 

18.25 

17.07 

16.28 

1.97 

0.92 

0.76 

0.56 

19.761 

13233 - 

1201 C 

3198 

16.18212 

2.05673 

370 

19.23 

17.68 

16.44 

2.79 

0.75 

0.70 

0.69 

19.675 

9057 C 


3144 

16.17702 

2.11280 

364 

18.61 

17.44 

16.40 

2.21 

0.88 

0.74 

0.48 

19.831 

12808 C 


1186 

16.22006 

2.15270 

310 

19.78 

18.31 

17.14 

2.64 

1.12 

1.04 

0.78 

20.448 

15883 - 


1093 

16.24090 

2.15469 

305 

19.50 

17.86 

16.56 

2.94 

2.23 

1.34 

0.92 

19.718 

15999 C 


3083* 

16.18528 

2.06219 

280 

19.05 

17.82 

16.79 

2.26 

0.71 

0.45 

0.54 

20.219 

9328 - 


3083* 

16.18528 

2.06219 

147 

19.05 

17.83 

16.96 

2.09 

0.64 

0.40 

0.86 

20.433 

9328 - 


4251 

16.23290 

2.10034 

263 

18.85 

17.54 

16.59 

2.26 

1.22 

1.07 

0.96 

20.040 

11806 C 



*There are two entries for G3083 with different periods as discussed in section 8. 


Table 7. Other large amplitude variables 


G 

RA 

equinox 

Dec 

2000 

J 

H K s 

(mag) 

A K S 

Sibbons+ 

Albert+ 

4098 

16.29459 

2.02230 

19.04 

17.60 

16.24 

1.0 

7773 C 


4165 

16.24627 

2.10189 

18.17 

16.98 

16.23 

0.4 

11937 C 

2671 C 

1183 

16.29165 

2.14605 

18.59 

17.51 

16.77 

0.4 

15435 C 

5994 C 

4051 

16.21875 

2.07353 

18.15 

17.03 

16.43 

0.8 

9922 C 

2173 C 

3100 

16.11923 

2.07968 

18.12 

17.11 

16.55 

0.5 

10310 C 

3047 C 

3068 

16.18036 

2.03580 

18.35 

17.34 

16.74 

0.4 

8218 C 

1037 C 

3247 

16.12183 

2.10549 

18.85 

17.93 

17.38 

0.43 

12218 C 

1892 C 

4200 

16.23416 

2.13410 

18.36 

17.25 

16.56 

0.7 

14536 C 

2843 C 

2111 

16.09110 

2.20471 

18.33 

17.14 

16.1 

0.40 

18795 C 




IC 1613 and its AGB Variables 13 


away from G1017 but which is almost certainly the same 
source. It is described as an LPV with V = 19.060, AV = 
1.87 mag (measurements in white light) and P=533.618 
days. 

G3011: Kurtev et al. (2001) found what they described 
as the first Mira variable in IC 1613, determined a period of 
640.7 days for it and assigned a spectral type of M3e. This 
is the source we call G3011 with the lithium-rich spectrum 
described in section 4. Note that the period of 640.7 days 
does not ht our JHKs or our R data We determined a 
period of 550 days, with an estimated error of about 30 days. 
It is not clear at this stage if the period has changed or if 
the earlier data could be consistent with a somewhat shorter 
period. It is also listed as an LPV, V2950, by Antonello et 
al. (2000) with a period of about 645 days and a white light 
amplitude of 1.7 mag. It will be important to establish if 
the period is changing as this could indicate that we are 
observing it at a time that is close to a helium shell flash. 

G1016: This was identified as V40, an irregular red 
variable, by Sandage (1971), and by Antonello et al. (1999), 
where it is V1872. It is listed as an Lc variable, 11613 V51, 
in the GCVS. 

On the basis of their variability these stars could be 
supergiants or AGB variables, but the spectrum discussed 
above (section 4) proves that G3011, and by implication the 
others, must be an AGB star undergoing hot bottom burn¬ 
ing. It seems likely that our survey is complete for these lu¬ 
minous O-rich Mira variables, even though it is very incom¬ 
plete for the faint M-type Miras and for the redder fainter 
C Miras. We would also not have detected highly reddened 
OH/IR sources, if such stars were present. 


9 COMPARISON WITH THEORETICAL 

ISOCHRONES 

Observations suggest that AGB stars end their AGB evo¬ 
lution as Miras, i.e. large-amplitude variables pulsating in 
the fundamental mode (Feast & Whitelock 1987, Whitelock 
& Feast 2000, Feast, Whitelock & Menzies 2002). Theory, 
(e.g., Marigo et al. 2008), suggests that there will be some 
evolution of the period and colour once fundamental pul¬ 
sation has started. However, the existence of the Mira pe¬ 
riod luminosity relation (Whitelock et al. 2008; Whitelock et 
al. 2013) together with kinematic evidence (e.g. Feast et al. 
2006) suggests that the evolution with period is very limited 
and/or very rapid at least among older populations, which 
also have low envelope masses. It is therefore our assump¬ 
tion that the Miras observed in IC 1613 are very close to the 
end of their AGB evolution. The situation with HBB stars 
is less clear and it is at least possible that they are overtone 
pulsators (Feast 2009), which will evolve into very long pe¬ 
riod, possibly OH/IR, stars when they become fundamental 
pulsators. 

Marigo et al. (2008) discuss the theoretical evolution 
of AGB stars with a wide range of ages, and they provide 
colours on the 2MASS system via a web interfac<[j Fig. [ 3 ] il¬ 
lustrates two representative isochrones from their work. The 
brighter of these passes close to the luminous M giants and 

1 stev.oapd.inaf.it/cgi-bin/cmd 


represents a relatively young population of just 2.5 x 10 8 yrs 
with Z=0.004. Stars on the upper part of this isochrone are 
indeed undergoing HBB. Increasing the metallicity moves 
the curve to the right, fitting lower luminosity M stars rather 
better, but missing the large amplitude variables. Decreas¬ 
ing the age moves it to higher luminosity. Assuming these 
M giant variables are the most luminous AGB stars then 
the AGB is observed down to 2 x 10 8 yrs. Note that for the 
purpose of illustrating these isochrones in Fig.[3]we omit the 
last few points where the colours vary dramatically over a 
very small range of initial mass. 

Marigo et al (2008) note that the bell-shaped profile in 
some of their theoretical isochrones (their hg. 1) is due to a 
combination of HBB and high mass loss. Their isochrone 
for z=0.004 and an age of 2.5 x 10 8 years clearly shows 
such a profile with the four luminous Miras lying near the 
brightness peak of the isochrone (see our Fig. 1). Their older 
isochrones do not show such a structure. This is consistent 
with the fact that we do not find any older (i.e. less luminous 
with shorter period) Miras showing evidence for HBB by be¬ 
ing above the PL relation. This is despite the fact that the 
galaxy contains such older stars. The Marigo et al. results 
for the isochrone in our Fig. 1 indicate that these four lu¬ 
minous Miras have an initial mass of 3.5 Mq and a current 
mass of 3.4Mq. They also indicate that slightly more ad¬ 
vanced stars on the isochrone are evolving rapidly with high 
mass-loss rates and are thus less likely to be found. These 
conclusions are not significantly changed if the metallicity is 
somewhat greater than z=004. This is shown by the z=0.008 
and z=0.019 isochrones in fig. 1 of Marigo et al. The initial 
(current) masses at the same age would be 3.6 (3.5) Mq for 
z=0.008 and 3.8(3.7) Mq for z=0.019. The presence of these 
four HBB stars indicates a significant population of age near 
2.5 x 10 s years. We might therefore expect a population of 
short period (1-2 day) Cepheids at this metallicity (see the 
Cepheid period-age relation of Bono et al. (2005)). Unfortu¬ 
nately such stars are at, or just below, the magnitude limit 
of the OGLE survey for Cepheids in IC1613 (Udalski et al. 
2001). However, a deeper survey of a small region of the 
galaxy (Bernard et al. 2010) shows a significant population 
of Cepheids with periods in the 1 to 2 day range. We con¬ 
clude that our HBB stars are restricted to small regions in 
the colour-magnitude and period-luminosity spaces because 
there is a lower limit to the age at which HBB can occur, 
coupled with the fact that there is only a sparse population 
of younger stars. 

The second, fainter isochrone represents a population 
with an age of 1 Gyr with Z=0.001 (the mean for IC 1613 
according to Skillman et al. 2014). It goes through the bulk 
of the AGB stars although it predicts a change from O- 
rich to C-rich at J — Ks < 1.0, bluer than the J — Ks ~ 
1.2 seen observationally. The pulsation periods predicted by 
the models do not bear any relationship to those actually 
measured. Nevertheless, the models have improved in recent 
years and seem to provide a plausible qualitative fit to many 
aspects of the data. 

The one low luminosity (presumed) O-rich star, G4237, 
with Ks = 17.3, must be significantly older than those 
discussed above and more like the Miras found in metal- 
rich Galactic globular clusters. If it is comparable to the 
AGB variables in 47 Tuc it is about 10 Gyr or possi¬ 
bly older (Garfa-Berro et al. 2014 and references therein). 
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Table 8. Red stars from Sibbons et al. not in the IRSF list; these 
are illustrated in Fig. [5] 


Sibbons+ 

RA 

Dec 

Jim 

M 

A '2 m 

8310 

16.203081 

2.037964 

20.581 

19.587 

18.368 

9409 

16.214640 

2.063861 

19.932 

18.087 

16.640 

10697 

16.304150 

2.085250 

19.618 

18.161 

16.531 

10922 

16.181959 

2.088694 

19.484 

17.841 

16.743 

11103 

16.091618 

2.091228 

19.025 

17.380 

16.291 

12772 

16.299379 

2.112383 

19.660 

18.030 

16.915 

13031 

16.219521 

2.115672 

19.007 

17.468 

16.088 

13166 

16.177704 

2.117122 

19.484 

17.792 

16.489 

13343 

16.252777 

2.119403 

19.201 

17.622 

16.759 

14597 

16.187973 

2.134992 

19.105 

17.662 

16.668 

16335 

16.108934 

2.159931 

17.915 

16.540 

15.543 

17381 

16.294380 

2.176472 

19.623 

18.531 

17.483 

18237 

16.172691 

2.192831 

19.464 

18.054 

17.015 

18859 

16.176949 

2.205914 

19.883 

18.500 

17.612 

18925 

16.248026 

2.207578 

19.240 

17.369 

16.524 

19795 

16.235384 

2.227403 

19.298 

17.805 

16.513 

20671 

16.235109 

2.252667 

20.377 

19.610 

18.194 

20772 

16.242823 

2.255681 

20.803 

19.089 

18.652 


The most luminous point on an isochrone of 10 Gyr and 
Z=0.004 (not illustrated) falls very close to this star in the 
colour-magnitude diagram. We anticipate that there would 
be significant numbers of such variables in IC 1613, with 
—7.5 < Mk < —6.5 (see Feast et al. 2002, fig. 3) and there¬ 
fore with 17.9 > Kg > 16.9. These stars will be fainter than 
our limiting magnitude during much of their pulsation cy¬ 
cle and it is not surprising that we do not actually identify 
others. 


10 BOLOMETRIC MAGNITUDES AND 
DISTANCES 

In the discussion of distances below we follow the same pro¬ 
cedure for determining bolometric magnitudes as in our pre¬ 
vious papers because this will give us consistent values that 
are useful, at least, for estimating distances via the period lu¬ 
minosity (PL) relation. Note that the caveats outlined at the 
start of section 7 of Whitelock et al. (2013) apply as much 
to the present work on IC 1613 as they did to NGC 6822. 

The photometry was reddening corrected assuming 
E(B — V) = 0.06 mag, which makes some allowance for 
reddening internal to IC 1613 (e.g. Georgiev et al. 1999). 

Bolometric magnitudes for presumed C stars were cal¬ 
culated by applying a colour-dependent bolometric correc¬ 
tion to the reddening-corrected K magnitudes on the SAAO 
photometric system (Whitelock et al. 2009). 

The Carpenter (2001; web page updat^]) formulation 
is used to transform from the 2MASS system (Table [6| to 
the SAAO system. The resulting bolometric magnitudes are 
listed in the same table. FigEH shows the PL relation. Nine 
stars is insufficient to determine both a slope and a zero- 
point for the PL relation, so we assume that the slope of the 


2 www.astro.caltech.edu/~j mc /2mass/v3/transformations/ 


PL-relation derived for the Large Magellanic Cloud (LMC) 
will apply to IC 1613 and use the expression; 

M bol = -3.31 log P + 3.90. (4) 

This enables us to determine the distance modulus of 
IC 1613 with respect to the LMC, assuming that the LMC 
has (m — M)o = 18.5 mag (Feast 2013). The fact that the 
slope of the PL relation determined for 50 stars in NCC 6822 
(Whitelock et a. 2013) was identical, within the uncertain¬ 
ties, to that of the LMC gives us confidence that this is a 
valid assumption. 

The distance modulus derived for IC 1613 is (m—M) o = 
24.37 ±0.08 mag, where the error is statistical and does not 
take into account any additional uncertainty in the slope of 
the PL relation or distance to the LMC. If we leave G3083 
out of the solution, because of the uncertainty in its period, 
we derive (m — M) o = 24.36 ± 0.09 mag. 

Whilst this is the first estimate of the IC 1613 distance 
from Mira variables, there have been many attempts to de¬ 
rive it in other ways, e.g., Classical Cepheids, RR Lyrae 
variables, TRGB, Red Clump. Summaries have been given 
by, among others, Bernard et al. (2010) and Scowcroft et 
al. (2013). Bernard et al. (their fig. 16) standardize previ¬ 
ous results to an adopted LMC modulus of (m — M) o = 
18.515 ± 0.085 and a reddening of E(B — V) — 0.025. 
They find a mean of various estimates to be (m — M)o = 
24.385 ± 0.014, where we have converted their value to our 
adopted LMC distance modulus (18.50 mag). The quoted 
error is internal only. Two recent multicolour studies of Clas¬ 
sical Cepheids in IC 1613 (Pietrzynski et al. 2006, Scowcroft 
et al. 2013), when adopting an LMC modulus of 18.50, lead 
to (m — M) o = 24.291 ± 0.035 and 24.31 ± 0.03 (internal 
errors only), respectively, and E(B — V) = 0.090 ± 0.019 
and 0.05 ± 0.01, respectively. Uncertainty in the reddening, 
illustrated by the three estimates given above, will obviously 
affect optical observations more severely than those in the 
near- or mid-infrared. Some modulus derivations, including 
those obtained from optical observations of RR Lyrae vari¬ 
ables, rely on an adopted metallicity. Others, including the 
Classical Cepheid determinations just mentioned, assume 
metallicity effects are negligible There is some evidence (Ma- 
jaess et al. 2014) that the Cepheids in the central region 
of IC 1613 are affected by crowding. Their results indicate 
that, allowing for such an effect by excluding Cepheids in 
the crowded centre, would increase the Cepheid distances 
just mentioned by ~ 0.06 mag. Finally, although a common 
LMC modulus is adopted, there is no guarantee that all 
the LMC calibrators are at the same mean distance (see for 
instance Feast et al. 2012). Given the above our current esti¬ 
mate of (m—M) o from Miras is well within the uncertainties 
of other estimates. As in our previous work on other local 
group galaxies, this confirms the reliability of these objects 
as distance indicators. 

The PL(iv) diagram is shown in Fig. |13[ where the line 
uses the distance modulus derived above in the relation for 
O-rich stars in the Galaxy (Whitelock et al. 2008): 

M k = —3.51 logP +1.10 (5) 

In comparison to the PL(A) relation in NGC 6822 
(Whitelock et al. 2013 fig 8) the absence of C stars below 
the PL relation is very noticeable. This is almost certainly 
a selection effect caused by the limited depth of our survey; 
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Log(P) 


Figure 12. Bolometric Period Luminosity Relation. O-rich and 
C-rich stars shown as open circles and closed circles respectively. 
The line is equation 4 with (m — M) = 24.37. 

many of the redder sources we found in NGC 6822 would not 
have been detected had they been at the distance of IC 1613 
and > 0.7 mag fainter. 

The fact that the luminous O-rich stars lie above the 
PL relationship in Figs |13| and |12| is another indication that 
these stars are undergoing HBB, supporting the conclusions 
from the presence of lithium (section 4) and the isochrone 
fits (section 9). The implications of this are discussed further 
in the conclusion (section 11) below. 


11 CONCLUSIONS 

It has already been established that star formation in 
IC 1613 took place over a very long period (Skillman et al. 
2014). Our observations support this and indicate a clear 
young (few xlO 8 yrs) AGB population, a significant popu¬ 
lation with ages up to a few x 10 9 yrs as well as indications 
of a much older (> 10 Gyr) population. A deeper survey is 
required to say more about the older AGB variables and to 
identify post-HBB stars with thick shells. 

It seems that IC 1613 offers a potentially very interest¬ 
ing environment in which to test AGB evolution and in par¬ 
ticular the lifetimes and flux contributions from HBB and 
stars with high mass-loss rates. Deeper and more detailed 
observations in the near-infrared and at longer wavelengths 
will be needed to provide an effective test of the models. This 
type of work could go some way to resolving the controversy 



Log(P) 


Figure 13. K Period Luminosity Relation. O-rich and C-rich 
stars shown as open circles and closed circles respectively. The 
line is equation 5 with (m — M) = 24.37. 


associated with the fraction of light from intermediate age 
AGB stars that is contributed to the integrated infrared flux 
of galaxies at moderate redshift (e.g., Maraston et al. 2006; 
Zibetti et al. 2013). 

Selecting stars by J — K$ colour rather than by their 
variability provides a distinctly different sample of AGB 
stars and in particular misses some of the most luminous 
large amplitude O-rich AGB stars. It will be important to 
recognize this in future. 

Whitelock et al. (2003) suggested that the luminous O- 
rich Miras that lie above the PL relation in the LMC are hot 
bottom burning. The discovery of lithium in the spectrum of 
G3011, a luminous Mira in IC 1613, together with the ages 
implied from the isochrone fits to this and similar variables, 
provides strong support for this idea. It also leads us to 
suggest that more attention should be given to stars with 
large amplitude pulsations that lie above the fundamental 
PL relation. These may be pulsating in the first overtone 
(Feast 2009) and, unlike fundamental pulsators, are not at 
the end of their AGB evolution. The Marigo et al. (2008) 
isochrones indicate that they have a great deal of mass to 
loose before becoming white dwarfs. It seem possible that 
they evolve into OH/IR stars, or possibly extreme C-stars, 
at a later stage after HBB terminates. As OH/IR stars they 
would have much longer periods, be enshrouded in dust, and 
presumably pulsate in the fundamental mode. 

This work demonstrates that a small (1.4m) telescope 
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can be used to make a very reasonable estimate of the dis¬ 
tances to Mira variables at around 750kpc. It will obviously 
be possible to use larger telescopes to get a bigger sample 
of AGB variables in IC 1613 and make a more accurate dis¬ 
tance measurement. 
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